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ABSTRACT 

New mid-infrared period-luminosity (PL) relations are presented for RR Lyrae variables in the 
globular cluster M4 (NGC 6121). Accurate photometry was obtained for 37 RR Lyrae variables using 
observations from the Infrared Array Camera onboard the Spitzer Space Telescope. The dispersion of 
M4’s PL relations is 0.056, and the uncertainty in the slope is 0.11 mag. Additionally, we established 
calibrated PL relations at 3.6 and 4.5 /im using published Hubble Space Telescope geometric parallaxes 
of five Galactic RR Lyrae stars. The resulting band-averaged distance modulus for M4 is jt = 11.399 ± 
0.007(stat) ± 0.080(syst) ± 0.015(cal) ± 0.020(ext). The systematic uncertainty will be greatly reduced 
when parallaxes of more stars become available from the GAIA mission. Optical and infrared period- 
color (PC) relations are also presented, and the lack of a MIR PC relation suggests that RR Lyrae 
stars are not affected by CO absorption in the 4.5 /im band. 

Subject headings: Stars: variables: RR Lyrae — Globular Clusters: individual: M4 — Stars: distances 
— Stars: horizontal branch 


1. INTRODUCTION 

RR Lyrae variables are important tools in the investiga¬ 
tion of many fundamental astrophysical problems. They 
provide crucial constraints on the physical mechanisms 
driving radial oscillations and their interplay with stellar 
evolution (Cox 1963; Christy 1966; Castor 1971). Fur¬ 
thermore, RR Lyrae stars offer the opportunity to study 
the morphology of the horizontal branch and the Oost- 
erhoff dichotomy (Oosterhoff 1939). Current empirical 
evidence indicates that the mean period of fundamental 
mode (RRab or FU) RR Lyrae stars in Galactic globular 
clusters (GGCs hereafter) shows a dichotomous distribu¬ 
tion at 0.55 (Ool) and 0.65 (OoII) days (Sandage 1993), 
where the latter group is more metal-poor. There is 
also evidence that the dichotomy is the aftermath of the 
hysteresis mechanism suggested by van Albada & Baker 
(1973), i.e. that the pulsation mode depends on the di¬ 
rection of the evolution inside the instability strip (Bono 
et al. 1997; Fiorentino et al. 2012; Kunder et al. 2013). 
Baade (1944) employed the RR Lyrae stars as a probe to 
identify the two main stellar populations in the Galaxy, 
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as well as to study the stellar content of the Galactic 
Bulge through low-reddening regions (Baade 1957; Plaut 
1968). 

Beyond stellar evolution, RR Lyrae variables have also 
played a key role in providing estimates of cosmological 
parameters, and have been instrumental in measuring 
the distances to a sizable sample of GGcs. This allowed 
the estimation of absolute cluster ages, and in turn set 
the lower limit on the age of the Universe (Hesser 1991; 
Sandage 1993; Buonanno et al. 1998; Marm-Franch et al. 
2009; Bono et al. 2010; VandenBerg et al. 2013). More¬ 
over, RR Lyrae variables have been used to estimate 
the primordial helium content using the A-parameter, 
i.e. the mass to luminosity relation of low-mass cen¬ 
tral helium-burning stars (Caputo et al. 1983; Sandquist 
2010 ). 

RR Lyrae stars are also the most commonly adopted 
Population II distance indicator. With a lower mass than 
classical Cepheids, they have the key advantage to be 
ubiquitous, and have been identified in both early and 
late type stellar systems (Van den Bergh 1999). Their 
individual distances can be evaluated using multiple di¬ 
agnostics, including a visual magnitude-metallicity rela¬ 
tion (Chaboyer et al. 1996; Bono et al. 2003; Cacciari & 
Clementini 2003) and a statistical parallax (Kollmeier 
et al. 2013; Dambis et al. 2013). More importantly, 
they obey to well defined near-infrared (NIR) period- 
luminosity (PL) relations (Longmore et al. 1986; Bono 
et al. 2001, 2003; Catelan et al. 2004; Braga et al. 2015). 
These PL relations extend to mid-infrared (MIR) bands 
where they have the potential of being very accurate dis¬ 
tance indicators due to lower extinction and smaller in¬ 
trinsic scatter (Madore et al. 2013; Klein et al. 2014). 
The use of the I,V-I reddening free period-Wesenheit 
(PW) relation to estimate the individual distances of 
RR Lyrae stars dates back to Soszynski et al. (2003) and 
to Majaess (2010). A more recent theoretical framework 
developed by Marconi et al. (2015) further supports the 
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use of optical, optical-NIR, and NIR period-Wesenheit- 
metallicity (PWZ) relations to determine individual dis¬ 
tances of RR Lyrse stars. Empirical validations to the 
above pulsation and evolutionary predictions have been 
provided by Braga et al. (2015) for RR Lyrae stars in the 
GGC M4 and by Coppola et al. (2015, in preparation) 
for RR Lyrae stars in the Carina dwarf spheroidal. 

The Carnegie RR Lyrae Program (CRRP) aims to take 
full advantage of the unique characteristics of these stars 
in order to reduce the remaining sources of uncertainty in 
the Hubble constant to ±2%. RR Lyrae MIR PL relations 
will be used as the foundation of an independent Popu¬ 
lation II cosmological distance scale to calibrate TRGB 
distances for nearby distances, which in turn can be used 
to calibrate Type la Supernova distances. To achieve 
this goal, we observed over 1,700 RR Lyrae variables in 
31 selected GGCs, as well as ~ 3,000 RR Lyrae stars in 
strategically distributed areas in the Galactic Halo and 
the Bulge. In addition, we observed 48 of the nearest, 
brightest and less reddened RR Lyrae stars intended to 
be used as zero point calibrators. These observations 
have been conducted during the warm mission of the 
Spitzer Space Telescope (Werner et al. 2004) Infrared Ar¬ 
ray Camera (IRAC, Fazio et al. 2004), at 3.6 and 4.5 fim 
wavelength, with a cadence designed to obtain complete 
coverage of each RR Lyrae variable over at least one full 
period. 

In this work we focus on the GGC Messier 4 (M4, 
NGC 6121). This cluster is an ideal laboratory for stel¬ 
lar population studies given its proximity to the Sun, 
which allows us to obtain accurate photometric and spec¬ 
troscopic data for member stars well below the main- 
sequence turnoff. Due to these characteristics, M4 has 
been the subject of intensive observational campaigns 
over a wide range of wavelengths. It has a well charac¬ 
terized differential extinction of E(B — V) = 0.37 ± 0.10 
mag, where the uncertainty is the dispersion due to dif¬ 
ferential reddening (Hendricks et al. 2012), and mean 
metallicity of [Fe/H] = —1.10, (Braga et al. 2015, and 
references therein). In support of the CRRP program, 
we have analyzed available multi-epoch optical and NIR 
data, leading to the identification and characterization of 
45 RR Lyrae variables (Stetson et al. 2014). From these 
observations, we have derived accurate distance moduli 
based on optical and NIR PL and PW relations (Braga et 
al. 2015). In this paper we combine our previous results 
with the new MIR data obtained as part of the CRRP 
campaign. 

In Section 2 we present our new Spitzer photometry. 
Light curves for all the RR Lyrae variables in our sample 
are measured in Section 3. In Section 4 we derive MIR 
PL and period-color (PC) relationships for the cluster 
RR Lyrae variables, while in Section 5 we calculate the 
M4 distance modulus by calibrating our PL zero point 
using five nearby calibrator RR Lyrae stars with known 
parallax, also observed as part of the CRRP program. 
Dependence of the PL zero point from metallicity is also 
discussed in Section 5, while Section 6 summarizes the 
results of this work. 

2. OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRY 

The results of our ground-based optical and NIR mon¬ 
itoring of the stellar population in M4 have been pub¬ 


lished in Stetson et al. (2014) and Braga et al. (2015). 
In this paper we extend our wavelength coverage to the 
MIR, by including multi-epoch 3.6 and 4.5 /jm photom¬ 
etry of a large portion of the cluster. 

Our IRAC observations of M4 were executed on 2013 
June 2 - 3 as part of the warm Spitzer Space Telescope 
Cycle 9 (PID 90088). The cluster was observed at IRAC 
3.6 and 4.5 /im bands at 12 equally spaced epochs over 
~ 14 hours. For each epoch we obtained 30 frames with 
a medium scale five-point Gaussian dither pattern and a 
3x2 mapping, taken using the IRAC full frame mode 
with a frame time of 30 seconds to provide a median 
total exposure of 804 seconds per pixel. An 11' x 9' area 
at the center of the cluster was covered in both the 3.6 
and 4.5 ^m bands. Due to the IRAC focal plane array 
configuration, in the outer regions of the cluster only one 
band is available, with the area just north of the center of 
the cluster only covered at 3.6 /im and the region south 
of the center only covered at 4.5 /im. The total area 
mapped by our observations covers a strip of 11' x 22', 
roughly cutting North-to-South through the center of the 
cluster, as shown in Figure 1. This is about one third of 
the total cluster area that was mapped in our ground- 
based optical and NIR observations, but is still ~ 1.5 x 
the size of the half-mass radius of M4 (3.'65 according to 
Harris 1996). 

The CRRP IRAC data were reduced starting from 
the Basic Calibrated Data (BCDs) generated by the 
Spitzer IRAC pipeline version SI9.1.0. Photometry 
was performed using the most recent version of the 
DAOPHOT/ALLSTAR/ALLFRAME suite of programs 
for crowded-field stellar photometry (Stetson 1987, 
1994). In order to mitigate the effect of blending in the 
photometry of the M4 crowded field, we used an input 
source catalog derived by the higher angular resolution 
optical and NIR images described in Stetson et al. (2014). 
In addition to higher resolution, this input source cata¬ 
log is also deeper. For the optical and NIR, accurate 
photometry was achieved down to 24 and 22 magnitudes 
respectively, compared to 17 mag for the MIR images. 
This catalog was matched to the Spitzer data after ap¬ 
plying higher order geometric corrections to the individ¬ 
ual IRAC BCDs than the default ones provided by the 
Spitzer IRAC pipeline. With this procedure we matched 
17,595 sources at 3.6 //m and 16,313 sources at 4.5 //m, 
for which we measured an instrumental magnitude with 
variable Point Spread Function (PSF) fitting. The astro- 
metric accuracy of our final source catalog is of the order 
of 0.3". 

For the optical and NIR images, good photometry was 
achieved down to 24 and 20 magnitudes respectively 

The photometry was calibrated to standard IRAC 
Vega magnitudes using aperture photometry on mosaic 
images. The mosaics were generated using the IRACproc 
package (Schuster et al. 2006), a Perl Data Language 
wrapper script for the standard Spitzer Science Center 
mosaicking software MOPEX (Makovoz et al. 2005). We 
adopted a mosaic pixel scale of 0.61"/pix, or half the 
native IRAC pixel scale. The stars we selected for aper¬ 
ture photometry calibration are far from the center of the 
cluster to reduce blending. In addition, we removed any 
possible variable stars, by avoiding sources with a Welch 
& Stetson variability index (Welch & Stetson 1993; Stet¬ 
son 1996) greater than ler from the median. We selected 
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a total of 38 stars at 3.6 /tin and 44 at 4.5 /im. We used an 
aperture of 3 native IRAC pixels (~ 3.66") with an annu¬ 
lus of 3 to 7 pixels (3.66" to 8.54"). The aperture correc¬ 
tions necessary to convert magnitudes derived with this 
aperture/annulus combination to the standard IRAC 10 
pixel aperture photometry are 1.1284 and 1.1274 for 3.6 
and 4.5 /m i respectively, calculated on isolated stars in 
the mosaics. The calibrating stars cover a range of 4 
magnitudes and are distributed randomly across the ar¬ 
ray, outside of the central region. 

Despite the relatively small size of our adopted pho¬ 
tometric apertures, the crowding of the field even out¬ 
side the center of the cluster (M4 is projected behind 
the Scorpius-Ophiuchus cloud complex) still results in 
blending contamination of aperture photometry for our 
calibration sources. These blended sources are however 
in many cases resolved by the PSF fitting procedure per¬ 
formed in DAOPHOT on the higher resolution visible 
and NIR images, from which our input catalog is de¬ 
rived. We thus used our point source photometric cata¬ 
log to improve the photometric zero point calibration by 
subtracting from the calibrator aperture photometry the 
flux (measured from the PSF-fitting photometry) of any 
blended sources found within our chosen aperture. For 
each blended source we estimated the contaminating flux 
by measuring the fraction of its flux falling within our 3 
pixel aperture, as a function of its centroid distance from 
the aperture’s center. For this task we used the IRAC 
PSF described in Hora et al. (2012), taking advantage of 
its higher sampling (~ 1/100 of the native IRAC pixels) 
to derive corrections for fractional pixel shifts. At 3.6 /im 
17 of our 38 calibrators were found to be contaminated 
by an average of 0.03 mag. The 4.5 /tm region was less 
crowded, resulting in only 7 out of 44 contaminated cal¬ 
ibrators, but contaminated by 0.11 mag on average. 

The magnitudes of the de-blended calibrators were 
then subtracted from their DAOPHOT instrumental 
magnitudes to derive the final calibration zero point. We 
found no trend in the calibration zero point with either 
brightness or array location. The uncertainty in our fi¬ 
nal calibration is 0.018 mag at 3.6 /iin and 0.017 mag 
at 4.5 /im, based on the standard deviation between the 
instrumental and measured magnitudes. This standard 
deviation of the calibration increases dramatically for the 
fainter stars in our calibration sample. If we select only 
calibrating stars brighter than 14 mag, then the error in 
the calibration is reduced to 0.015 at 3.6 /im and 0.013 
at 4.5 /tm. 

At the end of this process, we obtained a photometric 
catalog in both IRAC bands for all point sources detected 
in our field of view. From this catalog we matched 37 of 
the 45 RR Lyras stars with optical photometry derived 
by Stetson et al. (2014), 26 of them fundamental mode 
(FU, RRab) pulsators, and 11 pulsating in first overtone 
(FO, RRc). Of the remaining eight RR Lyrae stars in the 
cluster, seven are outside the area mapped in at least one 
of the two IRAC bands, and one (V64) is a field variable 
in the background of the cluster, too faint to be detected 
by IRAC. All M4 RR Lyrae stars are indicated in Figure 1 
by squares, with open red and filled blue symbols for FU 
and FO RR Lyrae stars detected in at least one IRAC 
band. 

The final magnitude of each non-variable star was de¬ 
rived by averaging the value measured in all epochs and 


in all dithers. For the RR Lyrae variables we instead pre¬ 
served the photometry for each epoch, averaging only on 
the individual dithers, after correcting the photometry 
in each dither for the warm Spitzer location-dependent 
photometric corrections available at the Spitzer Science 
Center web site. 10 The photometric uncertainty in each 
epoch is given by the standard deviation of all of the 
available measurements. A sample of the individual data 
points for one star, VI, is given in Table 1; all data 
is available in the online version of this article. A few 
RR Lyrae stars were covered in only one of the five dither 
positions, and one magnitude measurement was available 
in each epoch. For these stars (VI, V42, V2, and V29) 
the photometric uncertainty is instead derived from the 
repeatability parameter provided by DAOPHOT, which 
tends to be larger than the photometric error based ex¬ 
clusively on photon and background noise. 

3. RR Lyrae LIGHT CURVES AND AVERAGE MAGNITUDES 

Our MIR sample of M4 RR Lyrae variables consists of 
37 stars, whose location is shown in Figure 1. Of those, 
26 are pulsating in fundamental mode (FU, RRab) and 
11 are first overtone pulsators (FO, RRc). Not all sources 
have complete coverage in the IRAC dataset: of the 37 
RR Lyrae stars in the area mapped by at least one IRAC 
band, only 28 have photometry at both 3.6 and 4.5 /im, 
with six sources covered only at 3.6 /im and three covered 
only at 4.5 /mi. The periods of our RR Lyrae stars range 
between 0.2275 days (V49) to 0.6240 days (V39), plus a 
long-period FU RR Lyrae (V52) with a period of 0.8555 
days. 

Figure 2 shows representative color-magnitude dia¬ 
grams (CMDs) for M4 based on our MIR and optical 
photometry, with the IRAC 3.6 /tm magnitude vs. the 
B— [3.6] (top left ) and 7— [3.6] (bottom left) colors. Note 
how the high photometric quality of our data allows a 
well characterized turn-off in all color combinations, and 
a well populated sub-giant, RGB and AGB branch. The 
FO and FU RR Lyrae stars are represented by filled blue 
and open red squares respectively. We see a clear separa¬ 
tion of the FO and FU RR Lyrae stars in the horizontal 
branch of the CMD, with the FO bluer than the FU vari¬ 
ables (right panels). 

We derived MIR light curves for each RR Lyrae in 
our catalog by phasing their multiple-epoch photome¬ 
try based on the epoch of maximum (To) and the period 
derived by Stetson et al. (2014) using the Lomb-Scargle 
method from the visible data. Sample light curves are 
shown in Figure 3 for FO variables and 4 for FU variables. 
The light curves are shown in order of increasing period. 
Following Madore et al. (2013), for each star we derived a 
smoothed light curve using a Gaussian local estimation 
(GLOESS) algorithm, a method where a second-order 
polynomial is fit locally and points are weighted by their 
Gaussian distance from an interpolation point. A full 
description of the GLOESS method is given in Persson 
et al. (2004). The mean magnitude for all stars was then 
determined by integrating the smoothed light curve flux 
intensity over one period, then converting the result back 
into a magnitude. The variability amplitude in magni¬ 
tude was also derived from the smoothed light curve, and 
the uncertainty is calculated using only the data points 
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that lie within the Gaussian window of the maximum 
and minimum of the smoothed curve. 

The uncertainty in the mean magnitude is calculated 
by the sum in quadrature of photometric uncertainty of 
each observation (y / £of /N) and the uncertainty in the 
fit {A/(M\/Y2)) where N is the number of observations, 
(7 1 is the photometric uncertainty, A is the peak-to-peak 
amplitude, and M is the number of uniformly spaced 
points (following Scowcroft et al. 2011). The value of M 
ranges from five to twelve; five (e.g. V41) when our ob¬ 
servations cover more than two periods and twelve when 
our observations cover one complete period (e.g. V5). 
Some light curves however were randomly sampled (e.g. 
V52, where the observations do not cover a full period), 
and the uncertainty in the fit is instead A/\/12N where 
N is the total number of observations. Table 2 gives the 
3.6 and 4.5 fj ,m IRAC mean magnitudes and amplitudes 
for the 37 RR Lyras stars in the IRAC field. 

Figure 5 shows the luminosity amplitude vs. period 
(Bailey diagram) for all RR Lyrae stars: filled blue 
squares for FO and open red squares for FU pulsators. 
Black crosses indicate FU RR Lyras stars identified in 
Stetson et al. (2014) as Blazhko variables (the Blazhko 
effect is characterized by a conspicuous modulation of 
amplitude and phase; Blazko 1907). Recent long-term 
monitoring of Blazhko RR Lyrae by the Kepler satel¬ 
lite has provided evidence of not only period-doubling 
(Szabo et al. 2010) but also the possible occurrence of 
additional periodicities on a time-scale longer than the 
typical Blazhko period (Kolenberg et al. 2011; Chadid et 
al. 2011). As expected the FO and FU RR Lyrae variables 
are well separated in both period and amplitude, and the 
amplitude of the FU pulsators show a clear trend of de¬ 
creasing with period. Figure 6 shows instead the ampli¬ 
tude ratio between optical bands and the IRAC 3.6 pm 
band. The amplitude ratio plot can be used as a diagnos¬ 
tic for blends of stars with different color temperature. 
Two FO RR Lyrae stars (V40 and Cl) show an unusu¬ 
ally high optical to MIR ratio, possibly indicative of the 
presence of a faint optically-bright (hot) companion or 
blend. Another FO RR Lyrae (V49) shows instead an 
unusually low amplitude ratio, possibly the consequence 
of having an infrared-bright (cool) companion (or blend). 
It is also interesting to note that several (four out of 
seven at 3.6 pm and four out of six at 4.5 pm) of the 
FU RR Lyrae stars that lie below the general amplitude 
ratio trend are Blazhko variables, possibly due to these 
stars being observed at different phases in their ampli¬ 
tude modulation cycle (the optical and MIR data were 
not acquired simultaneously). If this is the case, the re¬ 
maining FU RR Lyrae stars in the same amplitude ratio 
group (V5, V3, and V52) could also be Blazhko stars 
that were missed by the period studies (or they could be 
stars with a faint cool companion or blend). 

4. MID-IR PERIOD-LUMINOSITY AND PERIOD-COLOR 
RELATIONSHIPS 

The left panels of Figure 7 show the IRAC PL relation 
for FO (filled blue squares) and FU (open red squares) 
RR Lyrae stars in M4. Blazhko variables are indicated by 
black crosses. The statistical uncertainties of the mean 
magnitudes are smaller than the size of the plotted sym¬ 
bols, and are not shown. Residuals from the PL relation 


were used as a diagnostic for blended variables. The 
central regions of M4 are not particularly dense, and no 
clear cutoff in radial distance from the center was appar¬ 
ent. Instead we employed a threshold in sigma to identify 
candidate blends. Two stars, V20 and V21, in the central 
regions of the cluster fall outside the threshold. The vari¬ 
able V20 was over 3cr above the fitted PL relation, and 
as shown in Figure 4, has a very poor quality light curve. 
The average magnitude of V21 was more than 5a above 
the fitted relation. These stars were not included in the 
final fit of the PL relations, and are shown in Figure 7 as 
open circles. 

All average magnitudes are corrected for extinction; 
albeit MIR extinction corrections are small compared to 
optical wavelengths, they still can not be ignored given 
the high photometric accuracy provided by IRAC. We 
used the total-to-selective extinctions of A[ 3 ^/E(B — 
V) = 0.203 and A [4 , 5] /E(B - V) = 0.156 from Monson 
et al. (2012) and a color excess of E(B — V) = 0.37±0.10 
from Hendricks et al. (2012). This results in extinction 
corrections of 0.075 ± 0.020 and 0.058 ± 0.016 mag at 3.6 
and 4.5 /irn, respectively. 

We calculated the zero point and slope (with their un¬ 
certainties and the best fit standard deviation) for the 
PL relation using an unweighted least squares fit. An 
unweighted fit is preferred to avoid biasing by brighter 
longer period variables, which tend to have smaller pho¬ 
tometric uncertainties. The PL relations take the form 

m = a + 6(log(P) + 0.55) (1) 

m = a + 6(log(P) + 0.26) (2) 

where log(P) = —0.55 and log(P) = —0.26 are represen¬ 
tative of the mean period of the FO and FU variables 
respectively. The scatter between the FO and FU PL 
relations is comparable. Due to the small number of 
FO variables, we also derived the global PL relations for 
all RR Lyrae stars (Figure 7, right), having fundamen- 
talized the periods of the FO variables using the relation 
log(PF[/) = log(PFo)+0.127 (Iben & Huchra 1971; Rood 
1973; Cox et al. 1983). The fundamentalized PL relation 
takes the form 

mn = a + 6(log(P) + 0.30) (3) 

where log(P) = —0.30 is representative of the mean fun¬ 
damentalized period of all RR Lyrae stars in our sample. 
The zero point a, slope b , their errors, and the standard 
deviations for all PL relations are given in Table 3. Note 
how the secondary modulation of the candidate Blazhko 
variables seem to have no measurable effect on their av¬ 
erage magnitudes, as they fit well on the PL relation. 

Figure 4 in Madore et al. (2013) showed that, us¬ 
ing available data, the slope of the RR Lyrae PL re¬ 
lation is a monotonically increasing function of wave¬ 
length, and asymptotically approaches a wavelength- 
independent slope. In Figure 8, we recreate this plot by 
adding the optical and near-infrared slopes from Braga 
et al. (2015), the IRAC MIR slopes calculated above, and 
the results from Klein et al. (2014). Our data is consis¬ 
tent with the values found by Madore et al. (2013) and 
Klein et al. (2014) for the WISE W1 and W2 bands, and 
confirm with a high level of accuracy that the slope ap¬ 
proaches the value of —2.60 at infrared wavelengths, as 
predicted by the period-radius relation in Burki & Mey- 
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lan (1986). 

4.1. Period-Color Relationship 

Figure 9 shows optical and infrared period-color (PC) 
relations for all M4 RR Lyras variables observed with 
IRAC (as before, filled blue and open red squares indi¬ 
cate fundamentalized FO and FU). We fit the data with 
an unweighted least squares procedure; the best fit pa¬ 
rameters are listed in Table 4. The figure clearly shows 
how the slope of the PC relation becomes more shallow 
as the distance between the two passbands decreases, as 
expected if the color is a temperature effect related to the 
radius of the pulsator. The V — K PC relation, in partic¬ 
ular, has been singled out for being the least affected by 
temperature uncertainties when deriving RR Lyras dis¬ 
tances with the Baade-Wesselink method (Carney et al. 
1992; Cacciari et al. 2000; McNamara & Barnes 2014), 
as opposed to using the standard B — V color. Our figure 
shows that PC relationships using MIR colors (and the 
IRAC 3.6 fim magnitude in particular) are even more 
advantageous, thanks to the steeper slope (i.e. ~ 1.93 
for the V — [3.6] relation vs. ~ 1.87 in the V — K fit) and 
slightly smaller dispersion around the best fit relation. 

The last panel in Figure 9 shows the IRAC [3.6] — [4.5] 
PC relation. In this case the slope is very shallow 
(0.046 ± 0.033), consistent with zero within ~ 1.4cr. The 
scatter around the best fit PC relation is also quite large 
(0.017 mag), three times larger than the uncertainty in 
the average magnitudes of individual stars. The lack 
of a significant [3.6] — [4.5] PC relation for RR Lyras 
stars suggests that these stars do not develop variable 
CO molecular absorption in their atmosphere, which is 
observed in other pulsating stars in the same instability 
strip (e.g. classical Cepheids). Scowcroft et al. (2011) 
found a well defined [3.6] — [4.5] PC relation for Large 
Magellanic Cloud Cepheids with logP < 2. The neg¬ 
ative slope of the PC relation and blue colors of these 
stars is explained by the dissociation and recombination 
of CO molecules in the IRAC 4.5 pan band, as previously 
noted by Marengo et al. (2010). This effect is particularly 
strong in Cepheids with P > 10 days, and is reflected in 
a precise phasing of the light and color curve of these 
stars. Shorter period Cepheids are instead characterized 
by a flat color curve and red colors. This is further sup¬ 
ported by Monson et al. (2012), which found that CO 
absorption drops off for temperatures greater than 6000 
K. Given that these RR Lyrae stars have higher effective 
temperatures, earlier spectral types and lower metallicity 
than Cepheids, it is not surprising that we do not find 
strong CO effects in the PC relation and color curve. 
The [3.6] — [4.5] PC relation appears to be flat, even 
when including Galactic RR Lyrae stars with very dif¬ 
ferent metallicities. Therefore, RR Lyrae observations at 

4.5 /im are not limited by strong metallicity effects as 
Cepheids are, and this band can be successfully used for 
distance scale measurements with similar accuracy as the 

3.6 /im band. The weakness of the 4.5 /im CO band in 
RR Lyrae stars, however, prevents this feature from being 
used as metallicity indicator for this class of variables. 

5. MIR PL ZERO POINT CALIBRATION AND M4 
DISTANCE MODULUS 

The CRRP project relies on establishing RR Lyrae vari¬ 
ables as highly accurate indicators for the first rung of 


the cosmological distance scale, by using their MIR PL 
relation. The zero point of this relation, however, needs 
to be calibrated. Efforts in this direction have been 
made in the recent past by using Galactic RR Lyrae 
stars observed by WISE (Madore et al. 2013; Klein et 
al. 2014), but the accuracy of the calibration has been 
limited by the uncertainties in the photometry and dis¬ 
tance of the calibrators. As mentioned in Section 1, the 
CRRP project addresses this issue with the observations 
of 43 Galactic RR Lyrae stars whose geometric paral¬ 
lax will be determined better than 2-3% with Gaia, and 
five bright RR Lyrae stars with parallax already deter¬ 
mined by Benedict et al. (2011) using the Hubble Space 
Telescope (HST) Fine Guidance Sensor (FGS). Ground- 
based monitoring and spectral observations of the cali¬ 
brators is ongoing, to ensure the characterization of their 
pulsation properties and, if required, to allow the calibra¬ 
tion of the metallicity dependence of the PL relation. 

The observation of a large sample of GGCs (themselves 
used as distance indicators) as part of the CRRP pro¬ 
gram provides a separate avenue to test the calibration of 
RR Lyrae stars and probe for metallicity effects. The dis¬ 
tance of M4, in particular, has been the subject of intense 
analysis, most recently by Braga et al. (2015) by using 
a theoretical and empirical calibration of the RR Lyrae 
optical and NIR PL relation for the cluster. The distance 
modulus derived in this work can be used as an alterna¬ 
tive calibration of the MIR PL zero point, as well as to 
test the effects of metallicity on the zero point derived 
from the five HST/FGS RR Lyrae calibrating stars. A 
similar approach will be followed for the remaining 30 
GCCs in the CRRP program (spanning a wide range in 
metallicity), for which independent distance moduli will 
be derived in the optical, NIR and MIR from their pop¬ 
ulation of RR Lyrae stars. 

5.1. Zero Point Calibration Using HST RR Lyrce Stars 

As mentioned above, at present there are only five 
RR Lyrae variables with available geometric parallax, 
i.e. the aforementioned stars observed with HST/FGS 
by Benedict et al. (2011): RR Lyr, UV Oct, SU Dra 
XZ Cyg (FU pulsators) and RZ Cep (FO). The four FU 
pulsators are the same stars used by Madore et al. (2013) 
to calibrate the RR Lyrae PL relation in the WISE bands. 
Even for these stars, however, the uncertainty in their 
distance modulus is still large (up to 0.25 mag), sig¬ 
nificantly reducing their effectiveness as calibrators for 
the distance scale. All five stars have been observed as 
part of the CRRP program in order to reproduce the 
Madore et al. (2013) zero point calibration using the 
IRAC 3.6 and 4.5 /im bandpasses, and with better phase 
sampling and much smaller (~ 10x better S/N ratio) 
photometric error. Their basic properties (derived from 
Benedict et al. 2011) are listed in Table 5: their peri¬ 
ods range from 0.3086 to 0.6642 days, and their [Fe/H] 
ranges from —1.43 to —1.83, significantly more metal- 
poor than M4 itself. A detailed analysis of the calibra¬ 
tion of the zero point using these five RR Lyrae stars 
at multi-wavelengths will be presented in a forthcoming 
paper. 

The mean magnitude of each star was computed as 
explained in Section 3. Each calibrator was then cor¬ 
rected for extinction, adopting the E(B — V) redden¬ 
ing from Benedict et al. (2011). The extinction correc- 
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tions were computed using the same A 3 , e /E(B — V) and 
A±. 3 /E(B — V ) relations as in Section 3. The absolute 
mean magnitudes, as well as the parallax, LKH correc¬ 
tion (Lutz & Kelker 1973; Hanson 1979), distance mod¬ 
uli, extinction corrections, and apparent mean magni¬ 
tudes are given in Table 5. Benedict et al. (2011) contains 
two conflicting results about the parallax to RZ Cep and 
although the value of 2.12 mas is their preferred solution 
(Benedict, private communication), we chose to adopt 
the 2.54 mas solution with an A v of 1 (assuming a mean 
spectral type of A5 (Preston 1959) and (B — V)o = 0.15) 
which appears to be more consistent with the new pho¬ 
tometric data. The dominating factor in the uncertainty 
of the absolute mean magnitudes is the uncertainty in 
their distances. 

Based on theoretical NIR PL-metallicity (PLZ) rela¬ 
tions, the slope of the RR Lyras PL relation is expected 
to be only weakly dependent on metallicity (Braga et 
al. 2015). Therefore, for this work we elected to adopt 
the well constrained slope of the fundamentalized PL re¬ 
lation derived in Section 4, rather than relying on just 
five Galactic calibrators (having a smaller range in pe¬ 
riod than the M4 RR Lyrae stars) to simultaneously fit 
both slope and zero point. Figure 10 shows the resulting 
PL relation, derived after fundamentalizing the period of 
RZ Cep. The calibrated PL relations derived using a one 
parameter fit are given by 

3 . 6] = —2.332(±0.106) log(P) - 1.176(±0.080) (4) 

a = 0.095 

M { 4 . 5] = —2.336(±0.105) log(P) - 1.199(±0.080) (5) 

er = 0.095 

where the uncertainty in the zero point is calculated by 
J]of/5 and eq are the uncertainties of the absolute 

magnitudes for the five calibrator stars used in the fit. 
The largest factor in the uncertainty of the zero point 
derived above is due to the uncertainty in the distance 
modulus of the calibrators, but it is a factor of 3 better 
than the value obtained in Madore et al. (2013) where 
both the slope and zero point were fit simultaneously and 
only four stars were used. In both IRAC bands, the slope 
we obtain is marginally shallower (but still consistent, 
within the error) than the values published by Klein et 
al. (2014) and Madore et al. (2013) for the WISE W1 and 
W2 bands. Further improvements in the uncertainty of 
the zero point will be possible when more accurate par¬ 
allaxes for all our calibrator RR Lyrae stars are obtained 
with the Gaia mission. 

5.2. M4 Distance Modulus 

From the zero point calibration of the MIR RR Lyras 
PL relation in the previous section, we can derive the 
true distance modulus of M4: 

jU [3 6 ] = 11.406 ± 0.010 (stat) ± 0.080 (syst) ( 6 ) 
±0.015 (cal) ± 0.020 (ext) 

^[4 5 ] = 11.391 ± 0.010 (stat) ± 0.080 (syst) (7) 
±0.013 (cal) ± 0.016 (ext) 


where we have adopted the M4 PL relation using the fun¬ 
damentalized period of FO pulsators. The error on the 
distance modulus comes from four sources. The statisti¬ 
cal error is the uncertainty in the zero point of the cluster 
PL relation, and is derived from the least squares fit to 
the IRAC data. The systematic error is the uncertainty 
in the calibrated PL zero point, given in the previous sec¬ 
tion. The calibration error is the dispersion in our pho¬ 
tometric zero point calibration to standard IRAC Vega 
magnitudes, described in Section 2. The extinction er¬ 
ror is the derived from the uncertainty due to differential 
reddening given in Section 4. 

The distance modulus we find is slightly larger than ex¬ 
pected, based on previous works (11.30±0.05 in Kaluzny 
et al. 2013, 11.28 ± 0.06 in Hendricks et al. 2012, and 
11.35 ± 0.03 ± 0.05 in Braga et al. 2015, this last value 
derived using RR Lyr itself as zero point calibrator), 
but agree within ltr. If we exclude the systematic er¬ 
ror, which will be significantly reduced once high accu¬ 
racy geometric parallaxes for a large sample of calibrators 
will become available from Gaia, we have an overall un¬ 
certainty better than 0.5% in each band. The average 
distance modulus of the two wavelengths then is: 

H = 11.399 ± 0.007 (stat) ± 0.080 (syst) ( 8 ) 
±0.015 (cal) ± 0.020 (ext) 

where the statistical error is reduced by a factor of ~ \/2 
(since by using both bands we take advantage of twice the 
number of available datapoints). Note that in this case 
the overall error is reduced to 1 %, small enough to allow 
the study of higher order parameters in the RR Lyras PL 
relation, first of all the possible dependence on metallic¬ 
ity, once data from more clusters is analyzed. 

Based on the models presented in Braga et al. (2015), 
while the dependence of the PL slope from metallicity 
is weak, we do expect a larger effect of [Fe/H] on the 
PL relation zero point. Therefore calibrating the PL 
relation with Galactic RR Lyrae stars, all with differ¬ 
ent metallicities, may induce a significant error. Fig¬ 
ure 11 addresses this issue, by plotting the residual be¬ 
tween the absolute mean magnitude of each calibrator 
and the value predicted by the PL relation, as a function 
of their metallicity. The metallicities used in the plot 
are listed in Table 5, and are the same listed in Bene¬ 
dict et al. ( 2011 ) (see references therein for the original 
source of the metallicity measurement), converted from 
the original (Zinn & West 1984) metallicity scale to the 
more recent Carretta et al. (2009) scale. For RR Lyr and 
UV Oct an extra step was necessary to correct for the dif¬ 
ferent solar abundances used by Kolenberg et al. (2010) 
when measuring the metallicity of these two stars, with 
respect to the Carretta et al. (2009). Figure 11 also plots 
the difference between the M4 distance modulus calcu¬ 
lated above for M4, and the value of 11.35 ± 0.03 ± 0.05 
derived by Braga et al. (2015) from NIR PL relations 
calibrated using RR Lyr. Due to the large uncertainties 
in absolute magnitude and metallicity of the Galactic 
RR Lyrae stars, no trend is apparent (slope of -0.2±0.3 
mag/dex). 

If we adopt the M4 distance modulus value of 11.283± 
0.001 ± 0.018 derived by Braga et al. (2015) from NIR 
data and theoretical PLZ relationships for RR Lyrae stars 
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(to avoid the large uncertainty in the calibrator distance 
that is affecting the zero point calibration), we obtain 
the following calibrated RR Lyrae PL relations for the 
two IRAC bands: 


Mj 3 . 6] = —2.332(±0.106) log(P) - 1.054(±0.020) (9) 
M[ 4 . 5] = —2.336(±0.105) log(P) - 1.091(±0.020)(10) 

An in-depth study of the MIR RR Lyrae PL relations 
based on theoretical models calculated for the IRAC 
bands will be presented in a forthcoming paper. 


6. CONCLUSIONS 

We have presented new RR Lyrae MIR PL relations for 
the nearby globular cluster M4. Accurate IRAC photom¬ 
etry allows us to reduce the error in the PL slope by a 
factor of at least two from previous works. We have also 
demonstrated that PL relations in the IRAC bands are 
consistent with previous calibrations using WISE pho¬ 
tometry. 

We also presented optical and infrared PC relations. 
The steeper slope and smaller dispersion of the V — [3.6] 
relation suggests that it could be preferred over V—K for 
use in the Baade-Wesselink method. Unlike for Cepheid 
variables, the [3.6] — [4.5] relation is very shallow, indicat¬ 
ing that there is no CO absorption, and therefore little 
metallicity effect, in the 4.5 (im band. As a result, the 
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4.5 £tm band can be used as a distance indicator. 

We calibrated the zero point of the MIR PL relations 
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The uncertainty in their distances is the largest source 
of error, and we will be able to provide a much more 
precise calibration using upcoming results from the Gaia 
mission. Due to the uncertainty of the available distances 
to Galactic RR Lyrae stars and the influence of metallic¬ 
ity in the PL relation, we provided a separate PL rela¬ 
tion, where the zero point is calibrated using a distance 
modulus to M4 based on theoretical NIR PLZ relations. 
Future work will include analysis the remaining GCCs in 
the CRRP program in order to provide a comprehensive 
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can then be compared to other distance indicators, such 
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Fig. 1.— Map of all M4 stars for which we have MIR photometric measurements. The open red and filled blue squares display the 
position of fundamental and first overtone RR Lyrae stars respectively. The black diamonds are RR Lyrae from our optical catalog that are 
not matched in the IRAC catalog. The solid cyan and magenta lines indicate the coverage of the IRAC 3.6 and 4.5 fim data, respectively. 
The dashed line indicates the approximate field of view of the optical photometry in Stetson et al. (2014). North is up, east is left. 
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Fig. 2.— Sample color-magnitude diagrams for all sources detected in the IRAC 3.6 band, with B and I photometry from Stetson et 
al. (2014). Fundamental and first overtone RR Lyrae variables are represented by open red and filled blue squares, respectively. The right 
panels show a close up of the horizontal branch. 
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Fig. 3.- IRAC light curves for all FO RR Lyrae variables in M4. The light curves are arranged from top to bottom in order of increasing 
period. The left and right panel display the light curves at 3.6 and 4.5 /im. The solid red line is the smoothed light curve generated using 
a Gaussian local estimator algorithm (see text for more detail). All light curves cover the same range in magnitude. 
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Fig. 4.— Same as Figure 3 but for all FU pulsators observed by IRAC. 
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Fig. 4 CONT.— Same as Figure 3 but for all FU pulsators observed by IRAC. 
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Fig. 5. — Plot of the IRAC band amplitudes versus period for R.R Lyras stars in M4 (Bailey diagram). First overtone and fundamental 
variables are represented by filled blue and open red squares. Candidate Blazhko stars are marked by black crosses. 




Fig. 6.— Ratio between the optical and MIR band amplitudes for the first overtone (blue filled squares) and fundamental (red open 
squares) variables as a function of period. The left plot shows the ratio of the B to [3.6] amplitude, while the right plot shows the ratio 
of V band to [3.6] amplitude. Black crosses identify candidate Blazhko variables. The two FO RR Lyrae stars (V40 and Cl) that show an 
unusually high optical to MIR ratio, could possibly indicate of the presence of a faint optically-bright (hot) companion or blend. The FO 
RR Lyrae (V49) with an unusually low amplitude ratio, could indicate an infrared-bright (cool) companion (or blend) 
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Fig. 7.— IRAC period-luminosity relations for M4 RR Lyrae stars. Left — FO (blue filled squares) and FU (red open squares) variables 
are fit independently (black solid lines). Right —- single fit with fundamentalized FO RR Lyras stars. All magnitudes are corrected 
for extinction. Random uncertainties in the average magnitudes of all variables are smaller than the symbols used, and are not shown. 
Candidate Blazhko stars are marked with black crosses. The open circles mark the positions of the two stars (V20 and V21) that were not 
included in the fit of the PL relation due to blends with nearby sources. 
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Fig. 8. — An updated version of Figure 4 from Madore et al. (2013) showing the increase of the slope of the RR Lyrae period-luminosity 
relation with wavelength. The optical and near-IR slopes are taken from Catelan et al. (2004), Sollima et al. (2006) and Braga et al. (2015). 
The MIR slopes from Klein et al. (2014), Madore et al. (2013) are marginally steeper than the values derived in this paper (using IRAC 
data only), but consistent within their uncertainties . For clarity, results from the same band are offset by 0.2 pm (Catelan +0.2, Sollima 
—0.2). Error bars are shown for all points, with the exception the theoretical prediction from Catelan et al. (2004). 
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Fig. 9.— The period-color relationship for the RR Lyrae stars in M4. FO and FU variables observed with IRAC are shown in filled 
blue and open red squares respectively, with fundamentalized periods for the FO RR Lyras. Optical and NIR photometry is derived from 
Stetson et al. (2014). The solid black line represents the unweighted least squares fit. Five Galactic RR Lyrae stars used for calibration are 
shown in the bottom panel as black circles (not included in the fit). All magnitudes are corrected for extinction. Note all panels except 
the bottom span the same range in magnitude. 
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Fig. 10.— The calibrated IRAC RR Lyrae period-luminosity relation. The filled circles represent the absolute magnitudes of five Galactic 
RR Lyrae stars. In order of increasing period: RZ Cep, XZ Cyg, UV Oct, RR Lyr (black point), and SU Dra. The period of RZ Cep has 
been fundamentalized with the same relationship adopted for the M4 FO RR Lyrae stars: log Pfu — log Pfo + 0.127. The PL relations 
are shown by the dashed lines, with a slope determined by the M4 cluster data. 
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Fig. 11. — Deviations from the period-luminosity relation as a function of metallicity. Galactic RR Lyrae stars are shown with solid circles 
(the darker point, with the smallest error bars is RR Lyr) and the average of all stars in M4 is shown as a open square. 
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TABLE 1 

Photometry for VI. 


MJD a 

[3.6] mag 

cr [3.6] 

MJD a 

[4.5] mag 

o-[4.5] 

56446.223 

11.354 

0.030 

56446.227 

11.267 

0.079 

56446.285 

11.376 

0.010 

56446.293 

11.317 

0.050 

56446.336 

11.394 

0.003 

56446.340 

11.318 

0.051 

56446.391 

11.329 

0.011 

56446.398 

11.261 

0.074 

56446.438 

11.302 

0.012 

56446.441 

11.224 

0.089 

56446.492 

11.312 

0.021 

56446.496 

11.248 

0.076 

56446.543 

11.358 

0.016 

56446.547 

11.285 

0.069 

56446.598 

11.385 

0.035 

56446.602 

11.312 

0.075 

56446.652 

11.344 

0.021 

56446.656 

11.268 

0.078 

56446.699 

11.309 

0.015 

56446.703 

11.240 

0.063 

56446.746 

11.295 

0.029 

56446.750 

11.233 

0.077 

56446.797 

11.340 

0.023 

56446.805 

11.278 

0.080 


Note. — Table 1 is published in its entirety in the electronic 
edition. A portion is shown here for guidance regarding its form 
and content. 


a MJD = JD - 2400000.5 days 
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TABLE 2 

Mean IRAC magnitudes and amplitudes for the RR LyrjE stars in M4. 


ID a 

a (J2000.0) 

5 (J2000.0) 

Period (days) 

log(P) 

[3.6] mag 

[4.5] mag 

Amp [3 . 6 ] 

Amp [ 4 5] 

Mode b 

VI 

16 23 13.67 
-26 30 53.6 

0.28888261 

-0.5392786 

11.278 ± 0.007 

11.244 ±0.021 c 

0.10 ±0.01 

0.09 ±0.04 

RRc 

V2 

16 23 16.30 
-26 34 46.6 

0.5356819 

-0.271093 

10.976 ± 0.027 c 

10.908 

± 

0.010 

0.25 ±0.06 

0.28 ±0.02 

RRab* 

V3 

16 23 19.47 
-26 39 57.0 

0.50667787 

-0.2952681 


10.982 

± 

0.008 


0.26 ±0.01 

RRab 

V5 

16 23 21.03 
-26 33 05.8 

0.62240112 

-0.2059296 

10.815 ± 0.012 

10.758 

± 

0.009 

0.15 ±0.04 

0.15 ±0.02 

RRab 

V6 

16 23 25.79 
-26 26 16.3 

0.3205151 

-0.4941515 

11.201 ± 0.013 




0.10 ±0.02 


RRc 

V7 

16 23 25.95 
-26 27 41.9 

0.49878722 

-0.3020847 

11.020 ±0.013 

10.977 

± 

0.010 

0.27 ±0.02 

0.28 ±0.01 

RRab 

V8 

16 23 26.16 
-26 29 41.6 

0.50822359 

-0.2939452 

10.941 ± 0.013 

10.896 

± 

0.009 

0.28 ±0.02 

0.27 ±0.01 

RRab 

V9 

16 23 26.80 
-26 29 48.0 

0.57189447 

-0.2426841 

10.869 ±0.012 

10.814 

± 

0.011 

0.26 ±0.02 

0.26 ±0.02 

RRab 

V10 

16 23 29.21 
-26 28 54.3 

0.49071753 

-0.3091684 

11.046 ± 0.015 

11.002 

± 

0.011 

0.32 ±0.03 

0.32 ±0.01 

RRab 

Vll 

16 23 29.98 
-26 36 24.5 

0.49320868 

-0.3069693 


11.029 

± 

0.010 


0.26 ±0.02 

RRab* 

V12 

16 23 30.82 
-26 34 58.9 

0.4461098 

-0.3505582 

11.160 ±0.023 

11.097 

± 

0.013 

0.34 ±0.04 

0.34 ±0.01 

RRab 

V14 

16 23 31.29 
-26 35 34.5 

0.46353111 

-0.3339211 

11.140 ±0.017 

11.083 

± 

0.014 

0.32 ±0.05 

0.31 ±0.04 

RRab* 

V15 

16 23 31.96 
-26 24 18.1 

0.44366077 

-0.352949 

11.170 ±0.014 




0.31 ±0.02 


RRab* 

V16 

16 23 32.50 
-26 30 23.0 

0.54254824 

-0.2655616 

10.880 ±0.011 

10.834 

± 

0.008 

0.24 ±0.03 

0.24 ±0.01 

RRab 

V18 

16 23 34.70 
-26 31 04.6 

0.47879201 

-0.3198531 

10.980 ±0.016 

10.896 

± 

0.014 

0.27 ±0.04 

0.25 ±0.03 

RRab 

V19 

16 23 35.05 
-26 25 36.3 

0.46781108 

-0.3299295 

11.101 ± 0.014 




0.27 ±0.03 


RRab 

V20 

16 23 35.39 
-26 32 35.9 

0.30941948 

-0.5094523 

10.953 ± 0.031 

10.901 

± 

0.031 

0.09 ± 0.06 

0.06 ±0.06 

RRc 

V21 

16 23 35.93 
-26 31 33.6 

0.47200742 

-0.3260512 

10.752 ± 0.016 

10.727 

± 

0.016 

0.22 ±0.04 

0.24 ±0.03 

RRab 

V22 

16 23 36.95 
-26 30 13.0 

0.60306358 

-0.2196369 

10.795 ± 0.014 

10.744 

± 

0.009 

0.18 ±0.04 

0.18 ±0.02 

RRab* 

V23 

16 23 37.33 
-26 31 56.1 

0.29861557 

-0.5248876 

11.053 ± 0.011 

11.043 

± 

0.013 

0.09 ± 0.02 

0.08 ± 0.04 

RRc 

V24 

16 23 38.04 
-26 30 41.8 

0.54678333 

-0.2621847 

10.922 ± 0.012 

10.900 

± 

0.010 

0.22 ±0.02 

0.23 ±0.03 

RRab* 

V25 

16 23 39.42 
-26 30 21.3 

0.61273479 

-0.2127275 

10.805 ± 0.013 

10.741 

± 

0.010 

0.26 ±0.04 

0.28 ±0.02 

RRab 

V26 

16 23 41.65 
-26 32 41.1 

0.54121739 

-0.2666283 

10.938 ± 0.013 

10.873 

± 

0.017 

0.32 ±0.02 

0.28 ±0.04 

RRab 

V27 

16 23 43.17 
-26 27 16.3 

0.61201829 

-0.2132356 

10.814 ±0.012 




0.24 ±0.03 


RRab 

V28 

16 23 53.60 
-26 30 5.30 

0.52234107 

-0.2820458 

10.984 ±0.013 

10.928 

± 

0.014 

0.36 ±0.02 

0.36 ±0.02 

RRab* 

V29 

16 23 58.25 
-26 21 35.1 

0.52248466 

-0.2819265 

10.977 ±0.013 c 




0.29 ±0.03 


RRab* 

V36 

16 23 19.45 
-26 35 49.0 

0.54130918 

-0.2665546 


10.900 

± 

0.011 


0.28 ±0.03 

RRab 

V37 

16 23 31.60 
-26 31 30.6 

0.24734353 

-0.6066994 

11.411 ± 0.013 

11.376 

± 

0.007 

0.06 ± 0.03 

0.05 ±0.01 

RRc 

V38 

16 23 32.87 
-26 33 03.5 

0.57784635 

-0.2381876 

10.755 ± 0.012 

10.708 

± 

0.011 

0.19 ±0.03 

0.20 ±0.03 

RRab* 

V39 

16 23 34.67 
-26 32 52.1 

0.623954 

-0.2048474 

10.823 ±0.012 

10.779 

± 

0.008 

0.17 ±0.03 

0.18 ±0.02 

RRab* 

V40 

16 23 34.59 
-26 30 44.2 

0.38533005 

-0.4141671 

10.875 ± 0.020 

10.822 

± 

0.011 

0.07 ±0.04 

0.08 ± 0.02 

RRc 

V41 

16 23 39.50 
-26 33 59.8 

0.2517418 

-0.5990447 

11.451 ± 0.011 

11.411 

± 

0.007 

0.09 ± 0.03 

0.08 ±0.01 

RRc 

V42 

16 24 02.00 
-26 22 13.0 

0.3068549 

-0.5130669 

11.286 ±0.032 c 




0.09 ± 0.06 


RRc 

V49 

16 23 45.25 
-26 31 28.4 

0.22754331 

-0.6429359 

11.470 ±0.011 

11.448 

± 

0.010 

0.05 ±0.03 

0.06 ± 0.02 

RRc 

V52 

16 23 24.06 
-26 30 27.8 

0.85549784 

-0.0677811 

10.488 ± 0.015 

10.428 

± 

0.015 

0.13 ±0.02 

0.14 ±0.02 

RRab 

V61 

16 23 29.76 
-26 29 50.3 

0.26528645 

-0.5762849 

11.433 ± 0.015 

11.361 

± 

0.006 

0.04 ± 0.03 

0.04 ±0.01 

RRc 

Cl 

16 23 36.43 
-26 30 42.8 

0.2862573 

-0.5432434 

11.244 ±0.029 

11.168 

± 

0.032 

0.07 ±0.06 

0.11 ±0.06 

RRc 


a The ID given in Clement’s catalog, with the exception of Cl, newly named in Stetson et al. (2014) 
k The pulsation mode (RRab=FU, RRc=FO). Candidate Blazhko stars are indicated by an asterisk. 
c The uncertainty is derived using the repeatability parameter given in DAOPHOT. 






22 


Neeley et al. 


TABLE 3 

Observed MIR period-luminosity relations for RR Lyr^: stars in M4. 


Band 

a a 

b a 

a 

a b 

b b 

<7 

a c 

b c 

(7 

[3.6] 

11.207 

±0.027 

FO 

-2.658 

±0.428 

0.079 

10.841 

±0.009 

FU 

-2.370 

±0.139 

0.040 

10.929 

±0.010 

FU±FO 

-2.332 

±0.106 

0.056 

[4.5] 

11.159 

±0.022 

-2.979 

±0.337 

0.057 

10.806 

±0.010 

-2.355 

±0.168 

0.045 

10.893 

±0.010 

-2.336 

±0.105 

0.054 


a PL parameters of the form m = a + 6(log(P) -f- 0.55) 
b PL parameters of the form m — a b(log(P) + 0.26) 
c PL parameters of the form m = a + &(log(P) + 0.30) 


TABLE 4 

Fundamentalized period-color relations. 



Color 


a a 



b a 


(7 

V 

- J 

0.808 

± 

0.018 

1.421 

± 

0.192 

0.102 

V 

- K 

0.985 

± 

0.021 

1.874 

± 

0.232 

0.123 

V 

- [3.6] 

1.021 

± 

0.021 

1.934 

± 

0.221 

0.117 

V 

- [4.5] 

1.047 

± 

0.023 

2.012 

± 

0.238 

0.121 

[3. 

6] - [4.5] 

0.034 

± 

0.003 

0.046 

± 

0.033 

0.017 


a PC parameters of the form color — a + fe(log(P) +0.3) 


TABLE 5 

IRAC MAGNITUDES FOR GALACTIC RR LyRJE STARS. 



RZ Cep 

XZ Cyg 

UV Oct 

RR Lyr 

SU Dra 

Period 

[Fe/H\ 

E(B - V) 

+.6] 

A [4.S] 

0.308645 
-1.80 ±0.2 
0.252 

0.051 

0.039 

0.466579 
-1.43 ± 0.2 
0.100 

0.020 

0.016 

0.542600263 
-1.56 ±0.11 
0.090 

0.018 

0.014 

0.566805 
-1.50 ±0.13 
0.042 

0.009 

0.007 

0.660419 
-1.83 ±0.2 
0.010 

0.002 

0.002 

Parallax (mas) 
LKH correction 
(m — M) o b 

2.54 ± 0.19 a 
-0.05 

8.03 ±0.16 

1.67 ±0.17 
-0.09 

8.98 ± 0.22 

1.71 ± 0.10 
-0.03 

8.87 ±0.13 

3.77 ±0.13 
-0.02 

7.14 ±0.07 

1.42 ±0.16 
-0.11 

9.35 ±0.24 

[3.6] 

[4.5] 

7.891 ± 0.007 
7.865 ± 0.006 

8.676 ±0.016 
8.645 ±0.017 

8.200 ±0.014 
8.174 ±0.014 

6.486 ± 0.014 
6.468 ± 0.014 

8.616 ±0.019 
8.588 ±0.019 

M[3.6] 

M\ 4.5] 

-0.190 ±0.160 
-0.204 ±0.160 

-0.324 ±0.221 
-0.361 ± 0.221 

-0.688 ± 0.131 
-0.710 ±0.131 

-0.663 ± 0.071 
-0.679 ± 0.071 

-0.736 ±0.241 
-0.764 ±0.241 


a Benedict et al. (2011) provides two different parallax values for this star. We adopted the value given in their 
Section 4.3.2, rather than the number listed in their Table 8, which appears to be inconsistent with the PL 
relation. 

k Our distance moduli are slightly different than the values in Table 8 of Benedict et al. (2011) due to typo¬ 
graphical errors in that paper. 














